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INTRODUCTION 


Important  aspects  of  this  analysis  are: 

a.  Two-parameter  distributions  are  dealt  with. 

b.  Both  parameters  are  estimated  from  the  sample 
(i.e.,  no  parameter  is  "assumed  to  be  known)*: 

c.  A  representative  sample  of  size  2  or  greater 
is  acceptable  data. 

The  results  of  the  analysis  will  be  presented  first 
for  the  sake  of  brevity. 

The  derivations  will  be  carried  out  for  the  normal 
population  only,  since  the  lognormal  case  is  only  trivi¬ 
ally  different. 

The  general  procedure  in  the  derivation  of  equation 
(1)  will  be: 


1. 

2. 

R  by 

3. 

f idence 

4. 


Find  the  jo 
Obtain  the 
making  a  chan 
Use  Neyman* 
interval  for 
Simplify  th 


int  densi 
density  an 
ge  of  vari 
s  method  t 

/f . 

e  resultin 


ty  of  and 

d  distribution  function  of 
ab  le  . 

o  obtain  a  one-sided  con- 

g  expressions. 
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NOTATION 


n  = 

ycc.cr- 
^  >1 

= 

X  = 

. 

A  - 

i»- 

>  - 
= 

c  = 

X^i/»  <<? 

>  *■ 

£  - 


sample  size 

normal  population  mean  and  standard  deviation 
estimators  Of and  ex' 
mission  life 
Population  reliability 


ft 


estimator  of 
4^1  -  standard  normal  dis- 

yr tt 

tribution  function 

standard  normal  fractile  function  (inverse  of 

distribution  function) 

estimates  of  sO 4,  and  & 

estimate  of  f\ 

confidence  level 

paramaters  of  the  lognormal  distribution 
estimates  of  (mean  and  standard 

deviation  of  the  logs  of  the  data) 

confidenced  reliability 
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3. 


RESULTS 


lity 


a.  Normal  population: 
The  equation  which 


yields 


confidenced  reliabi- 


The  inverse  problem  may  just  as  easily  be  solved, 
i.e.,  if  one  stipulates  desired  confidenced  reliability, 
one  can  solve  equation  (1)  for  the  necessary  mission 


life  x. 


b.  Lognormal  population: 

The  equation  which  yields  confidenced 


re liabi- 


(2) 
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The  foregoing  remarks  apply  equally  to  equation  (2) 

In  order  to  implement  the  numerical  solution  of 
these  equations,  the  following  programs  are  needed: 

1.  A  program  to  evaluate  the  standard  normal  dis¬ 
tribution  function . 

2.  An  integrating  program  (such  as  Romberg  inte¬ 
gration  )  . 

3.  A  univariate  nonlinear  equation  solver  (a  reli 
able  method  being  the  bisection  or  midpoint 


method ) . 

Needless  to  say , 
value  if  one  does  not 


thtse  equations  are  of  very  little 
have  a  digital  computer  available. 
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4.  DERIVATION  OF  RESULTS 

Let  ^  •  '  '  /v-7  )  be  normally  distributed 

with  mean an<*  standard  derivation  O"  , 

If  ^~PC  is  the  desired  mission  life,  population 

reliability  is  given  by: 


{?=  l-&( 


XX  - yCst 


a* 


(3) 


and  the  reliability  estimator  is  given  by 


y  * 

-r'  / 

-  l( 


y\ 

<T 


(4) 


/?  *  * 

Equation  (4)  gives  (/\  as  a  function  of  y*-^-*  and  o'  ; 

therefore  the  distribution  of  ^  may  be  obtained  if  one 

s\ 

first  finds  the  joint  distribution  of  y^~  and  &  • 

Form  a  bivariate  change  of  variable  using  the  dummy 
variable  ~2y  . 


*?-  /-i 


-K  - /M 


dr 


(4) 


=  O- 


(5) 
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Finding  the  inverse  of  this  relation: 


the  inverse  relation  is: 

^  _  ^  %,-A  «7) 

—  "ZX  (  8 ) 


now  , 


Where 


The  Jacobian  is  obtained 


oint  density 
first  f  using 


functions . 
equations 


(7)  and 


(8)  : 
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j-  _  9(A< >9)  __ 

9t,?>  k) 

9A  At 

9?  ^ 


oA 

9A 

2-a 

At 

9A 

9? 

piA 

9A 

99 


Sx  _  —  A 
"  aY%)^ 


( io) 


(ii) 
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but 


Therefore  from  equations  (9)  and  (13): 


(  12) 


(  13) 


10 


or 


t3 


/<£■*)**  *<?  ' 


(  14) 


After 


> 


*  *) 


is  determined,  may  b< 


integrated  out  (in  equation  (14))  to  obtain  the  density 
\ 

of  ft  . 

n°»,  is  normal  with  mean and  standard 

deviation  £T~  .  A  **. 

/-V7  dT” 

"M  =  — - -  is  chi-square  distributed 

'  '  J  <cr 

with  /?"7 —  /  degrees  of  freedom: 


The  density 


y  <? 


of  £T  is  therefore  given  by: 


^6 


z- 


but 


<T~ 


J?~*4  _  Z 
JSS-  cr~ 
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Hou  is  normally  distributed  with  mean  and 

'  A 

standard  deviation  ;  therefore,  the  density  of 

.  f^T 

is  given  by: 

L  (A_  —s* 

'2'7r 


cr 


/Ts\ 


(16) 


-1 


N  ow  , 
ab  les 


since  yCsL 
,  equations 


and 

(15) 


* 

and 


are  independent  random  vari- 
(16)  yield  the  joint  density  of 
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Substituting  this  expression  into  equation  (14): 


13 


Equation  (6) 


0— 


o~ 


(20) 


Therefore,  using  equations  (18),  (19),  and  (20): 


Now,  £J~  is  integrated  out  to  obtain  the  density  of  ^ 


14 


Since  ^  is  a  dummy  variable  of  integration  and  /?  is 
the  argument  of  ,  the  only  numbers  upon  which  the 

form  of  is  dependent  are  /n  and  e .  Therefore  , 

the  density  of  the  reliability  estimator  is  a  one  para¬ 
meter  (  ^  )  density  which  is  independent  of  the  failure 
density  population  parameters  and  mission  life. 

Changing  the  argument  of  h  to  avoid  confusion  and 
adding  the  subscript  $  to  h  to  indicate  its  depen¬ 
dence  on  the  population  reliability  ,  the  density  of 

^  is : 
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oO 

CZ 

-  /  7 


4  7^  <? 


(24) 


The  distribution  funct 


•  f  7  • 

ion  of  f\  l 


s  therefore  given  by: 


,  -  ^ 3-S} ± X-yJ^ 


'TlA'cJ' 

(25) 


The  Neyman  method  of  finding  a  one-sided  confidence 
interval  for  ^  may  be  explained  through  the  following 
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The  curve 


is  determined  by: 


The  information  in  the  failure  data  is  introduced  at  this 
point : 


then  determined  by  solving  equation 


(28) 


In 

of 


(28) 


order  to  solve 


<R 


given  by 


equation 
equat i on 


(28) 

(25) 


,  the  distribution  function 
should  be  simplified. 


r*\ 

sm  —  I  _  i. 


//(r)  -  j 

O  6 


(25) 
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Changing  the  order  of  integration: 


i-  I 


-  (g 
0 


But  from  equation  (12)  2_ 

-  r  2, 

?&r<2  y°< 

*  i  £.*■ 


18 


Substituting  (32)  and  (33)  into  (31) 


O 
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Separating  (34)  into  two  integrals: 


Le  t 


/"V» 


To  evaluate  I,  make  the  change  of  variable: 


=r 


—  ^ 


^  > 


—  /Wj 


(  36) 


£ 


(37) 


20 


but  from  ( 20  )  ,  K  = 


using  (38)  in  (35) 
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But  from  equation  (27), 


~PC_ - yCX 


Therefore  equation  (1)  is  obtained: 


Equation  (2)  is  obtained  analogously;  one  simply  thinks  of 
working  entirely  within  the  domain  of  the  logarithms  of 
the  data  points. 
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